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Abstract 

The tricritical behavior of the two-dimensional g-state Potts model with vacancies for 
< q < 4 is argued to be encoded in the fractal structure of the geometrical spin clus- 
ters of the pure model. The known close connection between the critical properties of the 
pure model and the tricritical properties of the diluted model is shown to be reflected in an 
intimate relation between Fortuin-Kasteleyn and geometrical clusters: The same transfor- 
mation mapping the two critical regimes onto each other also maps the two cluster types 
onto each other. The map conserves the central charge, so that both cluster types are in the 
same universality class. The geometrical picture is supported by a Monte Carlo simulation 
of the high-temperature representation of the Ising model (q = 2) in which closed graph 
configurations are generated by means of a Metropolis update algorithm involving single 
plaquettes. 



1 Introduction 



The two-dimensional g-state Potts models [1] can be equivalently formulated in 
terms of Fortuin-Kasteleyn (FK) clusters of like spins [2]. These FK clusters are 
obtained from the geometrical spin clusters, which consist of nearest neighbor sites 
with their spin variables in the same state, by laying bonds with a certain prob- 
ability between the nearest neighbors. The resulting FK, or bond clusters are in 
general smaller than the geometrical ones and also more loosely connected. The 
FK formulation of the Potts models can be thought of as a generalization of (un- 
correlated) bond percolation, which obtains in the limit q — > 1. The geometrical 
clusters themselves arise in the low-temperature representation of the pure model 
[3]. 

For q < 4, where the model undergoes a continuous phase transition, the FK clus- 
ters percolate at the critical temperature and their fractal structure encodes the corn- 
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Fig. 1. Trace possessing a double point. 

plete critical behavior. The thermal critical exponents are obtained using the cluster 
definitions from percolation theory [4]. With clusters and their fractal properties 
taking the central stage, the FK formulation provides a geometrical description of 
the Potts model. 

The concept of correlated bond percolation has been turned into a powerful Monte 
Carlo algorithm by Swendsen and Wang [5], and by Wolff [6], in which not individ- 
ual spins are updated, but entire FK clusters. The main advantage of the nonlocal 
cluster update over a local spin update, like Metropolis or heat bath, is that it dras- 
tically reduces the critical slowing down near the critical point. 

Although it was known from the relation with other statistical models that the phase 
transition of the Potts models changes from being continuous to first order at q = 4 
[7], initial renormalization group approaches failed to uncover the first-order na- 
ture for larger q. Only after the pure model was extended to include vacant sites, 
this feature was observed [8]. In a Kadanoff block- spin transformation, the vacant 
sites represent blocks without a majority of spins in a certain state, i.e., they rep- 
resent disordered blocks. In addition to the pure Potts critical behavior, the site 
diluted model also displays tricritical behavior, which was found to be intimately 
connected with the critical behavior [8]. With increasing q, the critical and tricrit- 
ical fixed points move together until at q = 4 they coalesce and the continuous 
phase transition turns into a first-order one. 

Recently, the cluster boundaries of two-dimensional critical systems have been in- 
tensively studied by means of a method dubbed "stochastic Loewner evolution" 
(SLE s )-a one-parameter family of random conformal maps, introduced by Schramm 
[9]. In this approach, the Brownian motion of a random walker is described by 
Loewner's ordinary differential equation, containing a random term whose strength 
is specified by a parameter R > [10]. Different values of R define different univer- 
sality classes. Various results previously conjectured on the basis of the Coulomb 
gas map [1 1,12,13] and conformal invariance [14] have been rigorously established 
by this method. The interrelation between SLE S traces and the Coulomb gas de- 
scription is made explicit by observing that the Coulomb gas coupling parameter g 
can be simply expressed by g — 1/R [15]. The nature of the SLE S traces changes 
with R: for < R < 1 the path is simple (non-intersecting), while for 1 < R < 2 
it possesses double points (see Fig. 1); and for R > 2 it is space-filling [16,17]. 
This change reflects a change in the critical systems described by the SLE S traces 
[15]: For 1 < R < 2 they represent the hulls of FK clusters of the g-state Potts 
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model with 4 > q > 0, while for \ < R < 1 they represent the closed graphs of 
the high-temperature representation of the 0(n) model with — 2 < n < 2 and at 
the same time also the external perimeters of FK clusters with dual parameter l/R 
[15,18]. For a recent overview from the mathematical point of view, see Ref. [18]. 

In this paper, the tricritical regime of the two-dimensional (annealed) site diluted 
g-state Potts model with < q < 4 is studied from the geometrical point of view. It 
will be argued that the tricritical behavior of these models is encoded in the geomet- 
rical clusters of the pure Potts model in the same way that the critical behavior is 
encoded in the FK clusters. The relation between geometrical clusters and tricriti- 
cal behavior was first established by Stella and Vanderzande [19,20] for the special 
case q = 2, i.e., for the Ising model. Using arguments based on renormalization 
group, conformal invariance, and numerical simulations, they showed that the geo- 
metrical cluster dimensions of the Ising model at criticality are determined by the 
q = 1 tricritical Potts model, as was earlier conjectured by Temesvari and Herenyi 
[21]. The values of two of the three leading tricritical exponents characterizing the 
geometrical clusters were already determined before by Coniglio and Klein [22]. 
The distinctive feature of the q = 1 tricritical model is that it is in the same univer- 
sality class as the Ising model defined by the central charge c = \. In addition, it 
has the same correlation length exponent v — 1 as the Ising model. The boundaries 
of geometrical Potts clusters were also already known to be in the same universality 
class as the tricritical model with the same central charge [23,13]. Formulated in 
terms of SLE S traces, the hulls of geometrical clusters are described at criticality by 
traces with \<R<\, thereby forming the geometrical counterpart of the FK hulls 
which are described by the traces with 1 < R < 2. It will be shown here that not 
just the boundary dimensions, but all the relevant fractal dimensions characterizing 
FK clusters are in one-to-one correspondence with those of the geometrical clus- 
ters, thus providing a physical picture for the close connection between the critical 
and tricritical behaviors just mentioned. The important aspect of the map is that it 
leaves the central charge unchanged, so that both cluster types and their boundaries 
are in the same universality class, characterized by the same central charge. 

To support this geometrical picture, we carry out a Monte Carlo study of the high- 
temperature representation of the 2-state Potts, or Ising model. We generate the 
high-temperature graphs by using a Metropolis update algorithm involving single 
plaquettes [24]. By duality, the high-temperature graphs, which are closed, form 
the hulls of geometrical spin clusters on the dual lattice. We thus simulate the geo- 
metrical hulls directly without first considering the clusters. From the geometrical 
properties of these graphs, such as their distribution, the size of the largest graph, 
and whether or not a graph spans the lattice, the fractal dimension of the hulls of 
geometrical clusters can be determined immediately. Our numerical result agrees 
with the analytic prediction by Duplantier and Saleur [25], which was derived using 
the Coulomb gas map. 

The paper is organized as follows. In Sec. 2, those aspects of the diluted g-state 
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Potts model are reviewed that are of relevance for the following, in particular its 
cluster properties. Section 3.1 discusses the various fractal dimensions characteriz- 
ing FK clusters. In Sec. 3.2, the results for FK clusters are transcribed to geometri- 
cal clusters, which in Sec. 4 are shown to encode the tricritical Potts behavior. The 
Monte Carlo results for the high-temperature representation of the Ising model are 
presented in Sec. 5, followed by a summary in Sec. 6. 



2 Diluted g-state Potts Model 



The diluted Potts model can be defined by the Hamiltonian [26]: 

-m = KYtfw - 1) + ■'EOW, - - # £(<W - 1), (i) 

(ij) (ij) i 

where f3 denotes the inverse temperature and the double sum Z)(ij> extends over 
nearest neighbors only. The first term at the right hand with coupling constant K is 
the pure g-state Potts model with spin variable a, = 1, 2, • • • , q at the ith site. On 
this a second Potts model with auxiliary spin variable Tj = 1, 2, • • • , s, coupling 
constant J, and ghost field H is superimposed. The reasons for this extension of 
the pure Potts model are twofold. First, in the limit s — > 1, the s-state Potts model 
describes (bond) percolation [2], which is naturally formulated in terms of clusters. 
The extension thus allows for the investigation of cluster properties of the original 
model when the limit s — > 1 is taken. Second, the diluted model has two fixed points 
for q < 4. In addition to the pure Potts critical point, it also has a tricritical point, 
depending on the value of the coupling constant J. Adding vacancies therefore 
leads to two distinct scaling regimes, each with its own critical exponents. Note that 
the limit s — > 1 is subtle as precisely for s — 1, the Hamiltonian (1) is independent 
of J and H, so that it reduces to the standard Potts model. 

The second term in the Hamiltonian (1) connects nearest neighbors of like spins 
(for which S auaj = 1) with bond probability 

V = 1 - e" J , (2) 

while the ghost field H in the last term acts as a chemical potential for the auxiliary 
spins in the first state, Tj = 1. A configuration of auxiliary bonds thus obtained 
can be represented by a graph on a restricted lattice, consisting of those sites of the 
original lattice that have at least one nearest neighbor with like spins (see Fig. 2). 

Following Fortuin and Kasteleyn [2], we can rewrite the partition function of the 
model (1) as [26] 



TOW, 

. (ij) 



£P 6 (1-P)" 6 I1 [l + (s - l)e- Hn *] , (3) 
{r} c(r) 
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Fig. 2. Left panel: Original lattice with the spin variables taking three different values 
(q = 3) indicated by different grey scales. Right panel: Reduced lattice consisting of sites 
having at least one nearest neighbor with like spin, i.e., of the same shading. Auxiliary 
bonds, taking two different values (s = 2), are laid with a certain probability within clus- 
ters of like spins. The resulting bond clusters represent sites which are all in the same spin 
(indicated by dots with the same shading) and auxiliary spin state (indicated by bonds with 
the same shading). The broken bonds within a spin cluster are indicated by broken lines. 

where {T} denotes the set of bond configurations specified by b bonds and b broken 
bonds between like spins, while C(T) denotes the clusters in a given bond config- 
uration T. Finally, n c is the number of sites contained in the cth cluster, where an 
isolated single site counts as a cluster, irrespective of whether it is on the restricted 
or the original lattice. This last observation follows from the absence of any ref- 
erence to the spin variable in the last term of the Hamiltonian (1). In the partition 
function, the limiting case s = 1, where the Hamiltonian (1) is independent of J 
and H, can be recovered by noting that, for a given spin configuration the sum 
S{r} £> b (l ~~ P) b °f tne probabilities of all possible bond configurations adds up to 
unity. 

As mentioned above, cluster properties can be extracted from the partition function 
(3) of the diluted model by taking the limit s — > 1. Specifically [26], 



1 dlnZ 



N ds 



^ b:v\i-vf n ^ Hn n = £4e-* n , (4) 

iV \{r} C(T) I n 



where A" denotes the total number of lattice sites and i n is the cluster distribution 
giving the average number density of clusters of n sites. The thermal average in- 
dicated by angle brackets in Eq. (4) is taken with respect to the pure g-state Potts 
model, i.e., the first factor in the Hamiltonian (1). The right hand is seen to be the 
generating function for clusters. By differentiating it with respect to the ghost field 
H higher momenta in the cluster sizes can be obtained. 

The pure Potts part of the theory is easily dealt with by noting that, for a given 
spin configuration, it gives a factor ( e A ) > where a denotes the number of nearest 
neighbor pairs of unlike spin. 

The two fixed points of the diluted model correspond to two specific choices of the 
coupling constant J [22,27] with K fixed at the critical temperature K c of the pure 
Potts model (see Fig. 3). The first choice is obtained by taking J = K. This case 
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Fig. 3. Phase diagram of the diluted model (1) in the K-p plane, with p = 1 — exp(— J) 
the bond probability. The arrows indicate the renormalization flow in the infrared. The FK 
fixed point governs the critical behavior described by the FK clusters, while the geometrical 
fixed point governs the tricritical behavior described by the geometrical clusters. 

is special because the factors e~ K arising from the first term in the Hamiltonian (1) 
can now be related to the bond probability (2) as 



-K 



= 1 



P, 



(5) 



and the partition function becomes 



z FK =£Ai-pr 6 n[i+c< 

{r} c(r) 



l)e 



-Hn c 



(6) 



where B is the total number of bonds on the lattice, B = b + b + a. The sum Z){<r} 
produced the factor q since each cluster can be in any of the q spin states. For s — 1, 
this partition function reduces to the celebrated Fortuin-Kasteleyn representation of 
the Potts model [2], 

(7) 



z fk = Y j p\i-p) 1 



{r} 



where iVc is the number of clusters, including isolated sites, contained in the bond 
configuration Y. The clusters seen in the limit s — > 1 encode the complete thermal 
critical behavior of the model and are frequently referred to as Fortuin-Kasteleyn 
(FK) clusters. For q — > 1, the partition function (7) describes standard, uncorrected 
percolation, where the FK clusters coincide with the usual percolation clusters. 



The second choice is obtained by taking the limit J — > oo where the bond proba- 
bility p tends to unity. The only clusters surviving this limit are those without any 
broken bonds between like spins (b = 0), which are the geometrical clusters. The 
partition function can be written in this limit as 



{r} 



-Ka 



n [i + - ^~ Hnc 

C(T) 



P c (q), 



(8) 



where the factor P c (q) is such that the product over the clusters Ilc(r) Pc(o) ='■ 
Pr(q) gives the number of different spin configurations for a given bond config- 
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uration T. That is, Pr(q) is the number of g-colorings of the geometrical clusters 
contained in T, where it is recalled that an isolated single site counts as a cluster. 
For s = 1, the partition function reduces to 

Z G = Y / e~ Ka M<l), (9) 
{r} 

which is nothing but the standard low-temperature representation of the pure Potts 
model [3]. For the Ising model (q = 2), each graph can be colored in two different 
ways, Pr(2) = 2, so that the coloring number becomes irrelevant. For uncorrelated 
percolation (q = 1), P r (l) = 1 and a = trivially, so that only one geometrical 
cluster remains, representing a fully occupied lattice [27]. 



3 Pure Potts Model 



3.1 FK Clusters 



Adapting similar notations [12], we parameterize the two-dimensional g-state Potts 
models as 

yfq = -2cos(tt/k), (10) 

with 2 > R > 1 so that the argument of the cosine takes values in the interval 
[n/2, it]. Special cases are: 

• tree percolation (q — 0, R — 2) 

• uncorrelated percolation (q — 1, R — |) 

• Ising model (q — 2, R — |) 

• g = 3, k = | 

• g = 4, R = 1. 

The parameter k is related to the central charge c via [14] 

while the correlation length exponent z/ and the Fisher exponent r] C are given by 
[11]: 

2 1 1 3_ 

u =oo — -' ^c = 2----/t. (12) 

3 2 — k k 4 

The latter determines the algebraic decay of the cluster correlation function Gc( x ) 
at the critical point: 

G c (x,x')~l/|x-x / | d - 2+ ^, (13) 

with d the number of space dimensions. Physically, Gc(x, x') gives the probability 
that sites x and x' belong to the same cluster. The subscript "C" is to distinguish the 
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thus defined exponent from the standard definition based on the spin-spin correla- 
tion function. The other exponents can be obtained from the two given in Eq. (12) 
using standard scaling relations. 

For ease of comparison with the more often used notation where the central charge 
is given in terms of a parameter m [14]: 

mym + 1) 

which in turn is related to q through 

Jq = 2 cos ( - \ = -2 cos ( m n) , (15) 

with 1 < m < oo, we note that the relation with R reads for Potts models: 

k = , m = - . (16) 

m k — 1 

Usually, only the first equation in Eq. (15) is given, we included the second to 
clearly see the relation with the R parameterization. 

The critical behavior of the Potts model is also encoded in the FK cluster distribu- 
tion i n given in Eq. (4), which near the critical point takes the form 

i n ~ n- Tc exp(-#n), (17) 

as in percolation theory [4]. The first factor, characterized by the exponent t c , is an 
entropy factor, measuring the number of ways of implementing a cluster of given 
size on the lattice. The second factor is a Boltzmann weight which suppresses large 
clusters as long as the parameter 9 is finite. When the critical temperature is ap- 
proached from above, it vanishes as 9 oc (T — Tc) 1 /' 70 , with a c a second exponent. 
The cluster distribution then becomes algebraic, meaning that clusters of all sizes 
are present. As in percolation theory [4], the values of the two exponents specifying 
the cluster distribution uniquely determine the critical exponents as 

<x = 2 , fa = , 7c = , 

<jq ac oc 

Vc = 2 + d -, v = — I » D c = 7, (18) 

t c - 1 da c t c - 1 

where the fractal dimension of the clusters is related to the Fisher exponent (12) via 

D c = \{d + 2-r 1 c). (19) 
Various exponents are given the subscript "C" to indicate that the cluster definition 
is used in defining them. For FK clusters, where in terms of R 

12k(2 - k) _ 3R 2 + 2AR + 4 

<7C ~ 3/€ 2 + 8/€ + 4' T ° ~ 3R 2 + 8R + 4' ( } 
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the cluster exponents coincide with the thermal ones. The cluster definition not 
always yields the thermal critical exponents. For example, when the critical behav- 
ior of a system allows for a description in terms of other geometrical objects such 
as closed particle worldlines or vortex loops, a related but different definition is 
required to obtain the thermal exponents from the loop distribution [28]. 

The various fractal dimensions characterizing FK clusters and the leading thermal 
eigenvalue y T = \jv read in terms of R [13,29,30] 



D c = 


1 + 


D u = 


1 + 


D EP = 


1 + 


D RB = 


1 - 


Vt = 


3- 



1 3 
1 — R 

2R 8 


(21a) 


K 

2 


(21b) 


1 

2R 


(21c) 


3 R 
2R + 2 


(21d) 


3 

r- 


(21e) 



Here, D c is the fractal dimension of the clusters themselves, and Du that of their 
hulls [31,32]. In the context of uncorrelated percolation, the hull of a cluster can 
be defined as a biased random walk [31]: Identify two endpoints on a given clus- 
ter. Starting at the lower endpoint, the walker first attempts to move to the nearest 
neighbor to its left. If that site is vacant the walker attempts to move straight ahead. 
If that site is also vacant, the walker attempts to move to its right. Finally, if also 
that site is vacant, the walker returns to the previous site, discards the direction it al- 
ready explored and investigates the (at most two) remaining directions in the same 
order. When turning left or right, the walker changes its orientation accordingly. 
The procedure is repeated iteratively until the upper endpoint is reached. To obtain 
the other half of the hull, the entire algorithm is repeated for a random walker that 
attempts to first move to its right instead of to its left. The hull of FK clusters is a 
self-intersecting path. As remarked in the Introduction, these hulls for 4 > q > 
correspond to the SLE S traces with 1 < R < 2. 



The fractal dimension D EP characterizes the external perimeter. Its operational def- 
inition [31] in the context of uncorrelated percolation is similar to that for the hull 
with the proviso that the random walker only visits (nearest neighbor) vacant sites 
around the hull of the cluster. From the resulting trace those sites not belonging to 
the perimeter, i.e., without an occupied site as nearest neighbor, are deleted (such 
sites can, for example, be visited by the walker when it makes a right or left turn on 
a square lattice). This leads to non-intersecting traces. 

These algorithms have recently been used in a numerical study carried out on fairly 
large lattices (L = 2 12 = 4096) to determine the fractal dimensions of FK clusters 
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of the g-state Potts models with q — 1, 2, 3, 4 [33]. 



By construction, the external perimeter is smoother than the hull, or at least as 
smooth, so that D EP < D H . The two are equal when the fractal dimension D RB 
of the so-called red bonds [34] is negative. (A red bond denotes a bond that upon 
cutting leads to a splitting of the cluster.) For FK clusters, D EP is strictly smaller 
than Dn for all 1 < R < 2, while for R = 1 (q = 4), the fractal dimension of 
the red bonds becomes zero, and the hull and the external perimeter have identical 
dimensions. 

The two boundary dimensions are seen to satisfy the relation [30] 

(L> EP -l)(L>H-l) = i (22) 

An analogous relation, this time also involving the fractal dimension of the FK 
clusters themselves reads 

(£ C -£ EP )(Z) C -§L> H -I) = A (23) 

In addition, the fractal dimensions satisfy the linear relation 

D C -D U = \(D EP -D RB ). (24) 

The dimension D c of the FK clusters approaches the number of available dimen- 
sions, D c — > 2 when q — > (R — > 2). In this limit, also the hulls of the FK clusters 
become space-filling, _D H — > 2. 

As in percolation theory [34], the fractal dimensions can be identified with renor- 
malization group eigenvalues of certain operators. For example, the fractal dimen- 
sion D c of the FK clusters coincides with the magnetic scaling exponent y H , 

D c = y H = d-(3c/v, (25) 

while that of the red bonds, -Drb, coincides with the eigenvalue y,j in the J direc- 
tion, 

£>rb = yj, (26) 

and therefore describes the crossover between FK and geometrical clusters [19,29]. 
Specifically, 1/y j determines the divergence of the correlation length £ at the crit- 
ical temperature K c when the bond probability (2) approaches the critical value 

p c = 1 - exp(- J c ), i.e., £ ~ (p c - p)' 1 ^-'. 

3.2 Geometrical Clusters 

Starting from the FK cluster dimensions (21), we next wish to obtain the analog ex- 
pressions for the geometrical clusters of the Potts model at criticality. As argued by 
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Vanderzande [27], both cluster types are characterized by the same central charge 
c. From Eq. (11) it follows that a given value of c does not uniquely determine R. 
Indeed, inverting that equation, we obtain two solutions for R: 

_ 13- C ±^/( C -25)(c-l) 

«± = 12 » < 27 > 

with k + > 1 and R_ < 1. The solutions satisfy the constraint 

= 1. (28) 

Since the right hand is independent of c, replacing R with leaves the central 
charge unchanged, c(R) — c(1/R). 

As an aside, note that the fractal dimensions of the hull and external perimeter of 
FK clusters are related by precisely this map [30], which is called duality in SLE 
studies. This duality, which in the Coulomb gas language corresponds to the earlier 
observed correspondance g — > 1/g (recall that g — 1/R) [15,18], is also at the root 
of the relation (22) between the two boundary dimensions. 

Applying the duality map R — > 1/R to the FK cluster dimensions listed in Eq. (21), 
we obtain 



D% = 1 + — (29b) 

£>rb = 1 + ^ - \* (29c) 

2/t=3-|t, (29d) 

where R = R + . These dimensions precisely match those conjectured by Van- 
derzande for geometrical clusters [27]. We therefore conclude that the geometrical 
clusters of the Potts model are images of the FK clusters under the map R — > 1/R 
for given R. Since this transformation at the same time maps the critical onto the 
tricritical regime (see below), it follows that the geometrical clusters (of the pure 
Potts model) describe the tricritical behavior in the same way as the FK clusters 
describe the critical behavior. For R — 1 (q — 4), the dimensions of the FK and 
geometrical clusters become degenerate, and the critical and tricritical behaviors 
merge. For q > 4, the phase transition is discontinuous [7]. The fractal dimen- 
sion (29a) was first given by Stella and Vanderzande [19] for the Ising case, and 
generalized to arbitrary 1 < R < 2 by Duplantier and Saleur [35]. 

The external perimeter dimension D EP of FK clusters has no image under the map 
R — > 1 /R. To understand this, recall that by smoothing the hull of a FK cluster 
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D c = l + i^ + \ (29a) 
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Table 1 

FK and geometrical fractal dimensions characterizing the g-state Potts models, with q = 



0,1,2,3,4. 



q 


c 


m 
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D c 
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4 
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15 
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5 
3 


11 

8 


1 


13 
24 


187 
96 


11 

8 


15 
8 


5 
8 


3 


4 

5 


5 


6 
5 


28 
15 


8 
5 


17 
12 


6 
5 


7 

20 


153 
80 


17 
12 


7 
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60 
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one obtains its external perimeter. For geometrical clusters, on the other hand, the 
fractal dimension of the red bonds is negative Drb < as follows from Eq. (29c) 
with 1 < R < 2, so that _D§ = D EP and the hull of a geometrical cluster is already 
non-intersecting. This is also reflected by the SLE« traces. Under the transforma- 
tion R — > 1/R, the self-intersecting traces with R > 1, representing the hulls of 
FK clusters are mapped onto simple traces with R < 1, representing the hulls of 
geometrical clusters [15]. 

For given R, the fractal dimension of the external perimeters of FK clusters coin- 
cides with that of the hulls of geometrical clusters, D EP = D^(= D^ p ). Since FK 
clusters are obtained from geometrical clusters by breaking bonds between nearest 
neighbors with like spins, the smoothing of the FK hulls apparently undoes this 
process again (as far as the boundaries are concerned). 

Table 1 summarizes the various fractal dimensions appearing in the Potts models. 
Note that for q = 0, the geometrical cluster dimension is larger than the number of 
available dimensions, D G > 2, making the geometrical clusters unphysical in this 
case. The equality D EP = D RB is typical for tree percolation. The value D G = 2 
for q = 1 agrees with the observation that in this case the geometrical cluster 
represents a fully occupied lattice [27]. The cluster dimension D C (R) possesses a 
minimum atl</t = 2/v / 3<|, allowing the models with q = 4 and q = 2 to 
have the same dimension. 

The physical meaning of the thermal eigenvalue y G is related to the existence of 
a critical magnetic field H S (K), such that in the region K < K c , H > H S (K) a 
geometrical cluster spanning the lattice (hence the subscript "s") is always present. 
The eigenvalue namely determines the vanishing of this field as K approaches 
K c from below [19]: 

H s (K)~(K c -KYt. (30) 
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4 Tricritical Potts Model 



In the traditional representation in terms of the parameter m determining the central 
charge c through Eq. (14), the g-state tricritical Potts models is parameterized by 
[14] 

Jq = 2 cos (— ) = -2 cos (^-^ttV (31) 

\mj \ m J 

where m is restricted to 1 < m < oo. As for the Potts models, usually only the 
first equation is given. The second one is included because on comparing with the 
parameterization (15) of the Potts models, we observe that the two are related by 
inverting m/(l + m). Given the connection (16) with the R notation, it follows that 
this is nothing but the central-charge conserving map R — > 1/R which relates the 
FK and geometrical clusters of a given Potts model. Whereas R + parameterizes the 
Potts branch (10), the solution R_ of Eq. (27) parameterizes the tricritical branch, 

y/q = — 2 cos(7t/k_), (32) 

with R_ restricted to the values | < k_ < 1, so that the argument of the cosine 
now takes values in the interval [it, 3tt] . (From now on, the R's are given a subscript 
plus or minus to indicate the solution larger or smaller than 1. Up to this point only 
the larger solution was used and no index was needed to distinguish.) The results 
obtained for the Potts models can be simply transcribed to the tricritical models, 
provided R + used on the Potts branch is replaced by /t_ . This close relation between 
the two models was first observed by Nienhuis et al. [8]. Note that for | < k_ < | 
the geometrical cluster dimension exceeds the available number of dimensions, 
which is unphysical. The eigenvalue y^f (29d) with R — R + = 1/R- is one of the 
two leading thermal eigenvalues of the tricritical Potts model [12]. The second one 
is given by the inverse correlation length exponent l/u, with 

v=--±—. (33) 
4 1 — K_ 



Because replacing R + with R_ is tantamount to replacing R + with 1/R+, it fol- 
lows that for given central charge, the FK clusters of the tricritical model are the 
geometrical clusters of the Potts model (and vice versa). For example, the fractal 
dimension Dq of the FK clusters in the critical regime, Dq = 1 + 1/ 2R + + 3R+/8, 
translates into 1 + l/2/c_ + 3R-/8 = 1 + R+/2 + 3/8R+ for the tricritical regime. 
This is precisely the fractal dimension (29a) of the geometrical Potts clusters. 

As q increases, the critical and tricritical points approach each other until they an- 
nihilate at the critical value q — 4, where the FK and geometrical clusters coincide 
and the red bond dimension vanishes. As stressed by Coniglio [29], the vanishing 
of D RB signals a drastic change in the fractal structure, anticipating a first-order 
phase transition. 
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5 High-Temperature Representation 



5. 1 Monte Carlo Study 

To support the picture discussed above we carry out a Monte Carlo simulation of the 
high-temperature (HT) representation of the 2-state Potts, i.e., Ising model, adopt- 
ing a new update algorith [24]. HT, or strong coupling expansions can be visualized 
by graphs on the lattice, with each occupied bond representing a certain contribu- 
tion to the partition function. For the Ising model, defined by the Hamiltonian 

-(3H = (3j2S i S j , S t = ±l, (34) 

(ij) 

where the coupling constant is taken to be unity, the HT representation reads [36]: 

Z=(cosh(3) 2N 2 N v\ (35) 
{r G } 

where {To} denotes the set of closed graphs specified by b occupied bonds, and 
v = tanh (3. Traditionally, HT expansions are carried out exactly up to a given or- 
der by enumerating all possible ways graphs up to that order can be drawn on the 
lattice. We instead generate possible graph configurations by means of a Metropolis 
update algorithm, involving single plaquettes [24] (see Fig. 4 for typical configu- 
rations generated in the low- and high-temperature phases). By taking plaquettes 
as building blocks, the resulting HT graphs are automatically closed-as required 
[36]. The update is such that all the bonds of a selected plaquette are changed, i.e., 
those that were occupied become unoccupied and vice versa [24] (see Fig. 5). A 
proposed update resulting in b' occupied bonds is accepted with probability 

p H T = min(iy'- 6 ), (36) 

where b denotes the number of occupied bonds before the update. With I denoting 
the number of bonds on the plaquette already occupied, b and b' are related through 

V = b + 4 - 21. (37) 

In the following, we focus exclusively on the graphs and measure typical cluster 
quantities, such as the graph distribution, the size of the largest graph, and whether 
or not a graph spans the lattice. From this, the temperature where the graphs pro- 
liferate as well as their fractal dimension can be extracted as in percolation theory. 
Both the proliferation temperature and the associated correlation length exponent 
turn out to coincide with their Ising counterparts. 

By the well-known Kramers -Wannier duality [37], the HT graphs form Peierls do- 
main walls [38] separating spin clusters of opposite orientation on the dual lattice. 
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□ if 

D 




Fig. 4. Typical graph configurations, somewhat resembling the oil painting Rhythm of 
a Russian Dance by De St ij I artist Theo van Doesburg (1883-1931), generated on a 
16 x 16 square lattice with periodic boundary conditions in the high- (left panel) and 
low-temperature (right panel) phase. 

selected plaquette 



Fig. 5. Update mechanism at work. Left panel: Present graph with the plaquette selected for 
updating indicated by the broken square. Right panel: New graph after the update proposal 
is accepted. Both the old and new graph consist of 8 bonds, in accordance with Eq. (37) 
since two bonds on the plaquette were already occupied. 



Fig. 6. A plaquette on the original lattice corresponds to a Peierls domain wall on the dual 
lattice, separating the site at its center with reversed spin from the outside. 

Each bond in a HT graph intersects a nearest neighbor pair on the dual lattice of 
unlike spins perpendicular to it. In other words, the HT graphs are the boundaries 
of geometrical spin clusters (albeit on the dual lattice) whose fractal dimension we 
wish to establish. The advantage of the plaquette update we use is that these bound- 
aries are simulated directly without first considering the corresponding cluster. At 
the critical temperature, the domain walls lose their line tension and proliferate. 



When interpreted as domain walls, the HT graphs should strictly speaking be cut 
at the vertices, so that the graphs break down in separate polygons without self- 
intersections that only touch at the corners where the vertices were located. How- 
ever, it is expected that this does not change the universal properties at criticality 
we wish to determine. 

From the duality argument it also follows that the plaquette update is equivalent 
to a single spin update on the dual lattice (see Fig. 6). To illustrate this, the inter- 
nal energy U is computed, using the plaquette update. On an infinite lattice, the 
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Fig. 7. Exact internal energy U on a 16 2 lattice for the Ising and dual model [transcribed to 
the original model using Eq. (38)], and the Monte Carlo data obtained using the plaquette 
update. 

Kramers-Wannier duality implies that observables calculated at an inverse temper- 
ature j3 in the original Ising model can be transcribed to those of the dual model 
at an inverse temperature (3. The relation between the two temperatures follows 
from noting that, an occupied HT bond represents a factor v, while a nearest neigh- 
bor pair on the dual lattice of unlike spin on each side of the HT bond carries a 
Boltzmann weight exp(— 2/3), so that [37] 

tanhp = e-^, (38) 

or sinh 2/5 = 1/ sinh 2/3. 

On a finite lattice with periodic boundary conditions, however, a mismatch arises 
because single HT graphs wrapping the (finite) lattice are not generated by the pla- 
quette update-such graphs always come in pairs. The HT Monte Carlo study will 
therefore not exactly simulate the Ising model with periodic boundary conditions, at 
least not for small lattice sizes. (For larger lattices, single graphs wrapping the lat- 
tice become highly unlikely, so that their absence will not be noticed anymore, and 
the HT Monte Carlo simulation becomes increasingly more accurate.) In contrast, 
on the dual side, where the plaquette update corresponds to a single spin update, 
this class of graphs is not compatible with the periodic boundary conditions, so that 
they should not be included. Hence, the plaquette update simulates the (transcribed) 
dual rather than the original model itself. Figure 7 gives the exact internal energy 
of the original Ising model on a finite lattice with periodic boundary conditions 
[39,40] and that of the dual model transcribed to the original one using Eq. (38). In 
the figure, also the data points obtained using the plaquette update are included and 
seen to indeed coincide with the dual curve. For increasing lattice sizes, the dual 
and the original curves approach each other. 
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Fig. 8. Probability Pg for the presence of a spanning graph as function of @ measured 
for different lattice sizes L. Within the achieved accuracy, the curves cross at the thermal 
critical point j3 = (3 C . 



5.2 Simulation 



To determine the graph proliferation temperature, the probability P s for the pres- 
ence of a graph spanning the lattice as function of (3 is measured for different lattice 
sizes [4]. For small f3, P s tends to zero, while for large j3 it tends to unity. We con- 
sider a graph spanning the lattice already when it does so in just one direction. 
Ideally, the curves obtained for different lattice sizes cross in a single point, mark- 
ing the proliferation temperature. It is seen from Fig. 8 that within the achieved 
accuracy, the measured curves cross at the thermal critical point, implying that the 
HT graphs (domain walls) lose their line tension and proliferate precisely at the 
Curie point. 

The data was collected in 3.3 x 10 5 Monte Carlo sweeps of the lattice close to 
the critical point and 1.1 x 10 5 outside the critical region, with about 10% of the 
sweeps used for equilibration. After each sweep, the resulting graph configuration 
was analyzed. Statistical errors were estimated by means of binning. 

Finite-size scaling [4] predicts that the raw Ps data obtained for different lattice 
sizes collapse onto a single curve when plotted as function of (/3//3 c — X)L 1 I V with 
the right choice of the exponent v. By duality, the relevant correlation length here 
is that of the Ising model, so that v takes the Ising value v = 1 . With this choice, 
a satisfying collapse of the data is achieved over the entire temperature range (see 
Fig. 9). 

Next, the cluster exponents o G and t g specifying the graph distribution, 

4 ~ b- TG e- e \ 9<x(f3- (3 c ) 1/aG (39) 
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Fig. 9. The raw data of Fig. 8 replotted as function of {(3/ [5 C — l)L l l u , with the Ising choice 
v = 1. The data collapse is satisfactory over the entire temperature range. 

Table 2 

Percolation strength and average graph size xc at the (inverse) critical temperature 
P c = ln(l + \/2)/2 = 0.440687- • • for various lattice sizes L (see Fig. 11). 



L 


16 


20 


24 


32 


40 


48 


64 


80 


96 


128 


-Poo 


0.1129(6) 


0.0986(7) 


0.0876(6) 


0.0744(7) 


0.0642(8) 


0.0569(8) 


0.0461(9) 


0.0426(11) 


0.0353(10) 


0.0312(12) 


XC 


8.90(4) 


10.44(4) 


11.91(5) 


14.86(8) 


17.44(10) 


20.03(16) 


24.55(22) 


29.61(34) 


33.13(39) 


41.93(81) 



are determined, where £ b denotes the average number density of graphs containing 
b bonds. To this end we measure the so-called percolation strength P^, giving the 
fraction of bonds in the largest graph, and as second independent observable the 
average graph size [4] (see Fig. 10) 



Xg 



(40) 



where the prime on the sum indicates that the largest graph in each measurement is 
omitted. Close to the proliferation temperature, these observables obey the finite- 
size scaling relations [41] 



P 00 = L-M"P{L/£), 



Xg 



L-»/*X(L/fl, 



(41) 



where £ is the correlation length and the critical exponents Pg^Ig are related to 
a G , t g through Eq. (18) written in terms of the variables appropriate for the graph 
exponents. Precisely at T c , these scaling relations imply an algebraic dependence 
on the system size L, allowing for a determination of the exponent ratios (see Table 
2 and Fig. 11). 

The data was fitted using the nonlinear least-squares Marquardt-Levenberg algo- 
rithm, giving 

Pg = 0.626(7) 7g = 0.740(4), (42) 
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Fig. 10. The percolation strength (top panel) and average graph size \G (bottom panel) 
as function of 0/f3 c f° r different lattice sizes L. 

with x 2 /d.o.f = 1.15 and 0.94, respectively, and where it was used that v 
These values are perfectly consistent with the fractions {3q — | = 0.625, 7q = 
0.75, leading to the exponents 



1. 



aG = iT 



27 

TG = TT' 



(43) 



and the fractal dimension 



D G 



11 



(44) 



of the HT graphs we were seeking. By duality, this fractal dimension equals that of 
the hull bounding the geometrical spin clusters (Peierls domain walls) on the dual 
lattice. Our numerical result agrees with Eq. (29b) with R = K Ising = | appropriate 
for the Ising model. 
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Fig. 11. Log-log plot of the percolation strength (top panel) and average graph size \G 
(bottom panel) at the critical temperature as function of the lattice size L. The straight lines 
0.642 L -0 - 626 and 1.139 L 740 , respectively are obtained through two-parameter fits. 



The value (44) was predicted by Duplantier and Saleur, using the Coulomb gas 
map [25]. Subsequent support for that prediction was provided by Vanderzande 
and Stella [20] who drew on earlier numerical work by Cambier and Nauenberg 
[42]. A first direct numerical determination was given by Dotsenko et al. [43]. Em- 
ploying the Swendsen-Wang cluster update, these authors analyzed the geometrical 
spin clusters and their hulls at the critical temperature. They extracted the fractal 
dimension from the resulting hull distribution, which is algebraic at the critical tem- 
perature. Apart from directly simulating the hulls with the plaquette update, another 
advantage of our approach is the use of finite-size scaling which is generally con- 
sidered more reliable than the extraction of exponents by fits to the data obtained 
for a fixed lattice size. 
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As argued in Sec. 4 for the general case, the geometrical Ising clusters correspond to 
the tricritical Potts model with R — R_ — 1/ Ri sing = |, which according to Eq. (32) 
is the tricritical q = 1 model [19]. In addition to the correlation length exponent 
v = 1 which we, in accordance with Eq. (33), observed numerically and the second 
thermal eigenvalue y^f = y, this tricritical behavior is further characterized by 
[12] /3c = <§j. This value follows from the scaling relation (25) with the fractal 
dimension Dq of the FK clusters replaced by that of their geometrical counterpart 
(29a) with k = «i sing = | [19]. 

The fact that the two correlation length exponents featuring in the critical q = 2 
and tricritical q — 1 Potts models are equal is special to this case, being a result 
of the Kramers -Wannier duality. Indeed, equating the correlation length exponent 
(12) of the critical Potts models and that of the tricritical Potts models given in (33) 
with k_ = 1/k to assure that both models have the same central charge, yields 
k = Rising = | as only physical solution. 



6 Summary 

In this paper, it is shown that the geometrical spin clusters of the pure g-state Potts 
model in two dimensions encode the tricritical behavior of the site diluted model. 
These clusters, formed by nearest neighbor sites of like spins, were shown to be 
mirror images of FK clusters, which in turn encode the critical behavior. Since the 
mirror map conserves the central charge, both cluster types (and thus both fixed 
points) are in the same universality class. The geometrical picture was supported 
by a Monte Carlo simulation of the high-temperature representation of the Ising 
model, corresponding to q — 2. The use of a plaquette update allowed us to di- 
rectly simulate the hulls of the geometrical clusters and to accurately determine 
their fractal dimension. 
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